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ABSTRACT 

If a source star is gravitationally microlensed by a multiple lens system, the result- 
ing light curve can have significant deviations fi'om the standard form of a single 
lens event. The chance to produce significant deviations becomes important when the 
separations between the component lenses are equivalent to the combined angular Ein- 
stein ring radius of the system. For multiple lens systems composed of more than two 
lenses, however, this condition is diSicult to meet because the orbits of such systems 
are unstable. Even if events are caused by a multiple lens system with stable orbits 
where a pair of lenses are closely located and the other component (third body) has 
a wide separation from the pair, photometrically identifying the lens multiplicity will 
be difficult because the event will be identified either by a binary lens event caused 
by the close pair lenses or a single lens event caused by the third body. In this paper, 
we show that if a seemingly binary lens event is astrometrically followed up by using 
future high precision interferometers, the existence of an additional third body can be 
identified via a repeating event. We show that the signatures of third bodies can be 
unambiguously identified from the characteristic distortions they make in the centroid 
shift trajectories. We also show that due to the long range astrometric effect of third 
bodies, the detection efficiency will be considerable even for third bodies with large 
separations from their close lens pairs. 
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1 INTRODUCTION 



If a source star is gravitationally microlensed by a multiple 
lens system, the resulting light curve can have significant 
deviations from the standard form of a single lens event. 
The chance to produce significant deviations induced by the 
lens multiplicity becomes important when the separations 
between the component lenses are comparable to the com- 
bined Einstein ring radius of the system, which is related to 
the physical parameters of the lens system by 
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where mtot ~ rrii is the total mass of the lens system, 
is the total number of the component lenses, rrii are the 
masses of the individual lenses, and Doi and Dos are the dis- 
tances to the lens system and the source star, respectively. 
Some fraction of binary lens systems can meet this require- 
ment, and dozens of candidate binary lens events have been 
reported by the lensing survey and foUowup observations 
(Udalski et al. 1994; Dominik & Hirshfeld 1994; Alard, Mao 
& Guibert 1995; Alcock et al. 2000; Afonso et al. 2000; Al- 
brow et al. 2000). 



However, for multiple lens systems that are composed 
of more than two lenses, it is difficult to photometrically 
identify the lens multiplicity. There are two reasons for this. 
First, multiple systems with all constituent lenses having 
separations between them equivalent to 6e are very rare be- 
cause the orbits of such systems are unstable. Second, even 
if events are caused by a multiple lens system with stable 
orbits where a pair of lenses are closely located and the other 
component (third body) has a wide separation from the pair, 
the close pair and the third body will behave as if they are 
two independent lens systems, causing the event to be iden- 
tified either by a binary lens or a single lens event. One case 
the lens multiplicity can be identified is when the source tra- 
jectory approaches closely both the close lens pair and the 
third body, causing a repeating event. However, due to the 
short range photometric elTect of the third body, the chance 
to produce repeating events is very low (see § 3). Although 
there was a claim that a triple lens system was discovered 
by Bennett et al. (1999), it was subsequently shown to be 
better explained by a rotating binary (Albrow et al. 2000). 
As a result, despite the substantial fraction of multiple sys- 
tems (Evans 1968; Batten 1973; Batten & Fletcher 1989; 
Fekel 1981; Mayor & Mazeh 1987), not a single multiple 
lens event has been reported to date. 
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Until now, lensing observations have been carried out 
only photometrically. However, by using several planned 
high precision optical interferometers, such as those to be 
mounted on space-based platforms, e.g. the Space Interfer- 
ometry Mission (SIM) and the Global Astrometric Interfer- 
ometer for Astrophysics (GAIA), and those to be mounted 
on very large ground-based telescopes, e.g. Keck and VLT, 
it will become possible to observe lensing events astromet- 
rically. When am event is astrometrically observed by using 
these instruments, one can measure the lensing induced dis- 
placement of the source star image centroid position with re- 
spect to its unlcnscd position (centroid shift S). Astrometric 
lensing observation is important because the lens mass can 
be better constrained with the measured centroid shift tra- 
jectory (Miyamoto & Yoshii 1995; Walker 1995; Paczyhski 
1998; Boden, Shao & van Buren 1998). 

In this paper, we show that astrometric lensing cam also 
be used for efficiently detecting third bodies in multiple lens 
systems. This is possible because due to the long range astro- 
metric effect of third bodies, the lens multiplicity of a large 
fraction of events can be identified via repeating events. We 
investigate the properties of the deviations induced by third 
bodies in the centroid shift trajectories of events and esti- 
mate the efficiency of third body detections expected from 
the future lensing observations. We note that this paper is 
extension of the work of Han ct al. (2002) who recently 
demonstrated the high efficiency of astrometric lensing ob- 
servations in detecting very wide binary companions. 



2 BASICS OF MULTIPLE LENSING 

If a source star is lensed by a multiple lens system, the lo- 
cations of the resulting images are obtained by solving the 
lens equation. When all lengths are normalized by the com- 
bined Einstein ring radius, the lens equation is expressed in 
complex notations by 



E 



mi/mtot 



(2) 



where Zi are the positions of the lenses, — £, + irj and 
z = X + iy arc the positions of the source and images, and z 
denotes the complex conjugate of z (Witt 1990). The mag- 
nifications of the individual images are given by the inverse 
of the determinant of the Jacobian of the lens equation eval- 
uated at the image position, i.e. 



|det J| 



det J : 



dz dz 



(3) 



Then the total magnification and the source star image cen- 
troid shift are obtained by 



Ni 



and 



8 = 



-c, 



(4) 



(5) 



where ^ — (5, 77) and Zj = {xj,yj) are the vector notations 
for the source and image positions and Ni represents the 
number of images. Since the lens equation describes a map- 
ping from the lens plane onto the source plane, to find the 



image positions for given positions of the source and the 
lenses, it is required to invert the lens equation. 

For a single lens system (A'^ = 1), the lens equation can 
be simply inverted. Algebraically solving the lens equation 
yields two solutions of the image positions, and the total 
magnification and the centroid shift are expressed in simple 
forms of 



A = 

and 
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where u is the dimensionless lens-source separation vector 
normalized by ^e- For a single lens event, the light curve 
has a smooth symmetric shape (Paczyriski 1986) and the 
centroid shift follows an elliptical trajectory (Walker 1995; 
Jeong, Han & Park 1997). 

For a multiple lens system (N > 2), on the other 
hand, the lens equation is non-linear and thus cannot be 
analytically inverted. However, since the lens equation can 
be expressed as a polynomial in z, the image positions are 
obtained by numerically solving the polynomial. For a N 
point-mass lens system, the lens equation is equivalent to a 
(A'^^ -|- l)-order polynomial in z and there are a maximum 
AT^ -|- 1 and a minimum AT -|- 1 images and the number of 
images changes by a multiple of two as the source crosses 
a caustic (Rhie 1997; Witt 1990). The caustic is the main 
new feature of the multiple lens system, which refers to the 
source position on which the magnification of a point source 
event becomes infinity. For a binary lens system {N = 2), the 
sets of caustics form close curves. For a multiple lens system, 
the caustic structure becomes more complex and can exhibit 
self-intersecting and nested shapes. The order of the poly- 
nomial rapidly increases as the number of lenses increases, 
and thus it becomes difficult to directly solve the equation 
for systems composed of many lenses. One commonly used 
numerical method that allows one to study the lensing be- 
haviors of multiple lens systems regardless of the number 
of lenses is the "inverse ray-shooting method" (Schneider & 
Weiss 1986; Kayser, Refsdal & Staboll 1986; Wambsganss 
1997). The disadvantage of using this method is that it re- 
quires large computation time to study the detailed struc- 
tures in the patterns of magnifications and centroid shifts. 

In our analysis, we investigate the lensing properties of 
triple lens systems instead of testing systems with various 
numbers of lenses. We note, however, that for systems com- 
posed of more than three lenses, the individual sets of lenses 
composed of the close pair lenses and each of the widely 
separated companion, in most cases, can be treated as in- 
dependent triple lens systems. In addition, for triple lens 
systems, one can directly solve the lens equation instead of 
using the very time-consuming inverse ray-shooting method. 



3 PROPERTIES OF MULTIPLE LENSING 
BEHAVIORS 

If a lens system contains an additional third body, the lens- 
ing behavior of the lens system is affected by the third body. 
The most important effect of the third body is that it makes 

the effective positions of other lens components shifted to- 
wards it. The approximate amount of the shift is 
gT r -i 
d Irsl' 



(8) 
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where rs is the position vector of the third body from 
the center of mass of the two close pair lenses (here- 
after binary center), d = |r3|/^E,B is the separation be- 
tween the third body and the binary center normalized by 
the combined Einstein ring of the two close pair lenses, 
^E.B = [4G(mi +m2)/c2][(l/D„i) - {1/D,,)]^^\ and gx = 
mz/{m\ + 7712) is the ratio of the third body mass to the 
total mass of the close pair lenses. If this effect is taken into 
consideration, however, the photometric lensing behavior of 
events during the time when the source passes the region 
around the close lens pair is well approximated by that of 
the binary lens events without the third body. Similarly, be- 
sides the inverse effective positional shift of the third body 
towards the close lens pair, the lensing behavior of events 
during the source's approach to the third body is well ap- 
proximated by that of the single lens event caused solely by 
the third body. That is, the close lens pair and the third 
body behave as if they are two independent lens systems. 
In addition, since the photometric effect of both the close 
lens pair and the third body are confined to narrow regions 
around the individual lens systems, even if an event is caused 
by a triple lens systems, the event, in most cases, will be 
identified either by a binary lens or a single lens event. 

However, if a seemingly binary lens event is astrometri- 
cally followed up, the existence of an additional third body 
can be identified with a significantly increased efficiency. 
This is because, compared to the photometric effect, the 
astrometric effect of the third body endures to a large dis- 
tance from it (Miralda-Escude 1996), and thus it can cause 
noticeable deviations in the trajectories of the source star 
image centroid shifts even when the source approaches the 
third body with a considerable separations. 

To show this, we construct maps of excess magnifica- 
tions and centroid shifts of an example triple lens system. 
The excess magnification and the centroid shift are defined 
respectively by 



_ At - Ab 
^~ At ' 

and 

A(5 = 5t — <5b. 



(9) 



(10) 



where At and St represent the magnification and the cen- 
troid shift of the exact triple lens system, while Ab and 5b 
represent those of the binary lens event without the third 
body. The constructed maps are presented in Figure 1. For 
the construction of the maps, we take the effect of the po- 
sitional shifts induced by the third body into consideration. 
The two close pair lenses of the triple lens system have a 
mass ratio of ge = m2/m\ = 1.0 and the separation be- 
tween them is ab = 1.0 in units of Ob.b- The third body has 
a mass ratio of qt = 0.6 and it is separated from the binary 

center by at ~ 15.0 (also in units of 6'e,b^ with an orien- 
tation angle of 45° with respect to the axis connecting the 
two close pair lenses. The locations of the lenses are marked 
by filled dots on the maps. Grey scales are used to represent 
the regions of deviations with e > 5%, 10%, and 20% for the 
excess magnification map and AS > 5%, 10%, and 20% of 



t According to the strategy of monitoring binary lens events for 
third body detections, the event caused by the close pair lenses 
will be the standard of all measurements. Therefore, we normalize 
all lengths and time scales in units o{ 9b,b B.nd Bi where b 
is the time required for the source to transit Se.b- 




Figure 1. The maps of excess magnification (upper panel) and 
centroid shift (lower panel) of an example multiple lens system. 
The lens system is composed of three lenses where a pair of 
lenses are closely located and the other component (third body) 
is widely separated from the pair. The two close pair lenses have 
a mass ratio of qb = m2/mi = 1.0 and they are separated by 
Ob = 1.0 in units of Se,b- The third body has a mass ratio of 
gx = rn3/{mi +m2) = 0.6 and separated from the center of mass 
of the close pair lenses by ax = 15.0 also in units of 6b,b- The 
filled dots represent the locations of the three lenses. Grey scales 
are used to represent the regions of deviations with e > 5%, 10%, 
and 20% for the excess magnification map and A6 > 5%, 10%, 
and 20% of 9^ b for the excess centroid shift map. The coordinates 
of the maps are centered at the center of the mass of the close pair 
lenses and all lengths are normalized by 9e,b- The straight lines 
with arrows represent the source trajectories of events whose light 
curves and centroid shifts are presented in Fig. 2. The solid curve 
near the close lens pair represents the caustics. The inset in each 
panel shows the enlarged view of the map in the region around 
the close pair lenses. The dashed circle in each inset represents 
the combined Einstein ring of the close pair lenses. 



^E,B for the excess centroid shift map. For a typical Galactic 
lensing event, the angular Einstein ring radius is ~ (0)10^ 
/i-arcsec, and thus AS = Q.19e,b corresponds to ~ (0)10 
/i-arcsec. For comparison, we note the SIM will be able to 
r neasure positional shifts of stars as small as several /x-arcsec 
(tittp:/ /sim.jpl.nasa.gov). 

From Fig. 1, one finds the following facts. First, due to 
the large separation between the third body and the close 
lens pair, both the photometric and astrometric effects of 
the third body are not important in the region around the 
close lens pair, implying that not only the photometric but 
also the astrometric lensing behaviors of triple lens events 
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Figure 2. Example centroid shift trajectories (left panels) and light curves (right panels) of events caused by the triple lens system 

whose maps of excess magnification and centroid shift arc presented in Fig. 1. The source trajectories responsible for the individual 
events arc marked in Fig. 1, where the number in each panel corresponds to that of the source trajectory. In each panel, the dotted and 
solid curves represent tlie centroid shift trajectories and light curves produced with and without the third body. The insets inside the 
left side panels show the detailed structures of the centroid shift distortions induced by the third body. The insets inside the right side 
panels show the enlargement of the part of the light curve during the source's passages around the region of the close lens pair. 



during the source's passage of the region around the close 
lens pair are well described by the binary lensing approxima- 
tion. Second, the region of significant astrometric deviations 
around the third body is much larger than the region of 
significant photometric deviations, implying that for an im- 
portant fraction of multiple lens events the existence of third 
bodies can be identified from extended astrometric followup 
lensing observations. 

To see the characteristics of the photometric and astro- 
metric deviations induced by the third body, in Figure 2, we 
present the light curves and the centroid shift trajectories of 
several events caused by the triple lens system. The source 
trajectories responsible for the individual events are marked 
in Fig. 1. In each panel, we also present the curves of the 
binary lens events computed without the third body taking 
the clfectivc positional shifts of the close binary lens pair into 
consideration (dotted curves). On the centroid shift trajecto- 
ries, we mark the centroid positions ('x' symbols) measured 
with a time interval of ^b.b /4 to show the changing rate of 



the centroid position. For an event with Je.b ~ 1 month, 
therefore, this interval corresponds to roughly a week. From 
the comparison of the centroid shift trajectories and the light 
curves of the corresponding events, one finds that the astro- 
metric signatures of the third body can be clearly identified 
from the characteristic loop-shaped distortions in the cen- 
troid shift trajectories, while the photometric signatures are 
too week to be noticed. One also finds that the astrometric 
deviations last for a long period of time. 



4 EFFICIENCY OF THIRD DETECTIONS 

In the previous section, we showed that astrometric lensing 
observations will enable one to identify the third body of a 
multiple lens system with little ambiguity from the charac- 
teristic deviations it makes in the centroid shift trajectories. 
In this section, we determine the efficiency of detecting third 
bodies expected from the future astrometric lensing obser- 
vations. For this determination, we use the formalism intro- 
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duced by Di Stefano & Scalzo (1999) and further developed 
by Han at al. (2002). 

If we define 6b and 63 as the smaUest separations from 
the source trajectory to the binary center and the third 
body, the orientation angle of the source trajectory with re- 
spect to the line connecting the binary center and the third 
body is represented by 



6b ± 63 



(11) 



where the sign '—'is for the case when the closest points on 
the source trajectory from the binary center and the third 
body are on the same side with respect to the line connect- 
ing the binary center and the third body, and the sign 
is when the closest points are on the opposite sides. If all 
lengths are normalized in units of 9e,b, equation (11) is ex- 
pressed by 



1 f P b ± y/qrPz 
d 



(12) 



where /3b and /Ss are the lensing impact parameters of 
the two independent events involved with the binary sys- 
tem composed of the close pair lenses and the single third 

body, respectively .0 Let us additionally define /Js^th such that 
among events that were identified to be affected by the close 
lens pair by approaching the binary center with an impact 
parameter /3b, only events with source trajectories passing 
the third body closer than /^s.th can be identified as multi- 
ple lens events from the centroid shift deviations larger than 
a threshold value A5th. Then, the fraction of multiple lens 
events whose third bodies can be astrometrically identified 
is computed by 



P = 



where the threshold orientation angles have values 

-1 ^-B + ^Wr 133, til ' 



a+,th = sm 
a_,th = sin" 



for d> + ^/gr/^s.th, 
I a+xh = vr/2 



Q_ th — sm 



-1 / i9b- yWfa.th \ 

V " ) 



for I /3b - y?T/33,th I < d < /3b + ^/qrPi 



and 



. a-,th = -7r/2, 



(13) 



(14) 



(15) 



(16) 



for d < I /3b — •JqrP'i.tiil- We note that the threshold orien- 
tation angles have difi^erent values depending on the relative 
size of the astrometrically effective region of the third body 
compared to the separation between the third body and the 
binary center. For a given value of the threshold centroid 
shift deviation, A5th, the corresponding threshold impact 
parameter to the third body is obtained by 

t By definition, the lensing impact parameter represents the clos- 
est separation between the source trajectory to the lens (the cen- 
ter of mass of the lens system) normalized by the (combined) 
angular Einstein ring radius of the involved lens (lens system), 
and thus /Jb = 6b/^^e,b and jSz = b3/9E,3, where 9b, 3 represents 
the Einstein ring radius of the third body. 




Figure 3. The expected efficiency of detecting third bodies from 
astrometric foUowup lensing observations as functions of the sep- 
aration, d, and the mass ratio, Qti between the third body and 
the close lens pair. We assume that a third body is detected if 
the deviation induced by the third body is larger than 10% of the 
combined angular Einstein ring radius of the monitored binary 
lens event, 9^ b- The separation is expressed in units of b- 



/33,th = 



2 VA5th/fe,3 



(A<5th/fe,3) 



-8 



(17) 



If one expresses the threshold astrometric deviation in terms 
of the fraction of the combined Einstein ring radius of the 
close pair lenses, / = A(5tii/^E,B, equation (17) becomes 



/3 



(18) 



In Figure 3, we present the determined efficiency in 
the parameter space of the separation and the mass ratio 
between the third body and the close lens pair, P(d, gx). 
For this computation, we set / = 0.1. For the impact pa- 
rameter to the close lens pair, we adopt /3b = 0.0, but we 
note that the dependency of the efficiency on the adopted 
value of /3b is not important as long as /3b ^ d. We also 
note that can be larger than 1.0 because the third body 
can be heavier than the total mass of the close pair lenses. 
From the figure, one finds that the efficiency is substan- 
tial even for third bodies that are widely separated from 
their close lens pairs. It is worth to note that the scal- 
ing between d and qt are quite linear. This can be under- 
stood as follows. The probability of detecting third bod- 
ies is proportional to the threshold orientation angle, i.e. 
P oc a oc siri^^ (y^gT/33^th/d) [see eq. (12) and (13)]. In the 
limit of d > ^g7/33,th, sin"^(^g7/33,th/d) ~ ^gT/33,th/d, 
and thus P oc y/qrPs^th/d. Since /33,th oc ^/qr [see eq. (18)], 
one finds that P oc qT/d. Therefore, the scaling between d 
and qT is linear. 
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5 CONCLUSION 

In this paper, we show that future astrometric lensing ob- 
servations will provide an efficient method of detecting third 
bodies of multiple lens systems, which could not have been 
detected from conventional photometric lensing observa- 
tions. We showed that the deviation in the centroid shift tra- 
jectory induced by the third body of a multiple lens system 
has a characteristic loop that can be clearly distinguished 
from other types of deviations and thus can be unambigu- 
ously identified. In addition, since the deviations last for a 
long period of time, detecting third bodies will be possible 
even from sparse astrometric sampling. 

This work was supported by a grant (BSRI-01-2) from 
the Basic Science Research Institute of Chungbuk National 
University. 
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